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We investigate the extrinsic spin Hall effect in the electron gas model due to transition- 
metal impurities based on the single-impurity Anderson model with orbital degrees of free- 
dom. Both the skew scattering and side jump mechanisms are analyzed in a unified way, 
and the significant role of orbital degrees of freedom are clarified. The obtained spin Hall 
conductivities are in proportion to the spin-orbit polarization at the Fermi level (Z • s)^ as 
is the case with the intrinsic spin Hall effect: skew scattering term cr|1j oc {I ■ s)fi,5\(jxx, and 
side jump term (jgjj oc {I ■ s)^, where 5\ is the phase shift for p {I = 1) partial wave. Further- 
more, the present study indicates the existence of a nontrivial close relationship between the 
intrinsic term crg^ and the extrinsic side jump term cTgjj. 

§1. Introduction 

There is growing interest in the spin HaU effect (SHE), which is the phenomenon 
that an electric field induces a spin current in a transverse direction. After its theo- 
retical prediction by Hirsch in 1999, which was originally proposed by D'yakonov 
and Perel' in 1971,^^ investigations were initially focused mostly in semiconductors. '^^ 
Recently, however, experiments have made great progress on electrical spin injection 
and detection in metals. Therefore, SHE in various metals has received consid- 
erable attention owing to its fundamental interest as well as its potential application 
in spintronics. 

The intrinsic SHE arises from the Berry phase of the multiband Bloch func- 
tion.^)'®) In d-electron systems, the spin Hall conductivity (SHC) takes a large value 
due to the "orbital Aharonov-Bohm phase factor" induced by the d-angular momen- 
tum with the aid of the atomic spin-orbit interaction (SOI).^-*"^^) In Ref. 12), the 
present authors pointed out, for the first time, that the intrinsic SHC is in propor- 
tion to the spin-orbit polarization at the Fermi level (Z • s)^. In accordance with the 
sign change of {I ■ s)^, SHC changes its sign. The sign of the spin Hall angle (SHA 

Q^SH = % -T^ ) calculated in this study is consistent with the recent experimental ob- 

servation based on the spin pumping method. Furthermore, more recently, SHC 
experimentally measured by Morota et al.^^^ is semiquantitatively consistent with 
the result in Ref. 11). This fact strongly indicates that the intrinsic mechanism 
is dominant for the experimentally observed SHE in 4d and 5d transition metals. 
In /-electron systems, moreover, which possess larger angular momentum, a giant 
intrinsic contribution arises from the same mechanism as d-electron systems.^^^'^''-' 

In addition to the intrinsic mechanism, SHE is also caused by impurity scattering 
with the aid of the SOI, which is known as extrinsic SHE.^)'^)'^®)"^^) A controversial 
study on the extrinsic SHE in Au/FePt was reported in Ref. 22). To understand this 
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experimental observation, Guo et al. studied the skew scattering mechanism due to 
Fe impurities in Au.^^^ Therein, they considered that the giant SHE is originated 
from the orbital-Kondo effect of Fe impurities in Au host. However, Mihajlovic et al. 
reported the smah SHA in Au for a 60 nm-thick Au Hah bar.^^^ Consistently with 
the latter reported, {I ■ s)^ of Fe will be small since the orbital magnetic moment is 
almost quenched because of the small SOI of Fe atom.^^-*"^^^ Independently of Ref. 
23), the present authors studied the skew scattering term due to the /-electrons. 
Therein, we found that the giant SHC is realized only in the case of a large spin-orbit 
polarization. 

More recently, Seki et al?'^^ and Sugai et a/.^^-* have studied the Au layer 
thickness dependence of the SHA, and concluded that the giant SHA originates from 
the the scattering on the surface. Later, Gu et al. studied the possibility for a surface- 
assisted skew scattering on Pt impurities. In another group, in addition, Niimi 
and Fert et al. have recently reported a large SHA in Cu due to Ir impurity. ^^^'^^^ 

In the extrinsic mechanism, not only the skew scattering mechanism but the 
side jump mechanism also exists. Compared with the skew scattering mechanism, 
the side jump mechanism is less well understood. To elucidate the major role of the 
orbital degrees of freedom in the extrinsic SHE, studies on both the skew scattering 
and side jump terms based on the single-impurity Anderson model (SIAM) are highly 
required. Moreover, the side jump term originates from the anomalous velocity, as 
is the case with the intrinsic term. This fact could indicate the existence of a close 
relationship between the side jump term and the intrinsic term.^^^ 

In this paper, we study the extrinsic SHE based on the SIAM for transition 
metal atoms. The analytical expressions for both the skew scattering term Ugy and 
side jump term are derived. As is the case with the intrinsic term (Tg|^, both (t||j 
and (Tgjj are proportional to the spin-orbit polarization at the Fermi level {I ■ s)^: 
•^SH ^ ' ^)iJ.^i^xx, and fjgjj oc {I ■ s)^, where 6i is the phase shift for p {I = 1) partial 
wave. Therefore, (Z • s)^ ~ 0(1) will be a necessary condition for a giant extrinsic 
SHE, which is realized in the rare-earth atoms. In contrast, {I ■ s)^ is small in 4d 
and 5d transition metals, which makes the extrinsic SHE small. Both the intrinsic 
and extrinsic terms can be discussed from the same viewpoint that is the spin-orbit 
polarization at the Fermi level (Z • s)^. 

In the present study, it is found that the contribution of the side jump term due 
to the impurity such as Pt in Cu host is comparable to that of the intrinsic term in 
Pt. In contrast, the condition for the skew scattering term to be huge is difficult to 
be realized, which not only {I ■ s)^ but also 5i should be large. 

§2. Model and Hamiltonian 

In this study, we use the SIAM for transition metal atoms. In the presence of 
the SOI, d-electron states are specified by the total angular momentum J = 2 it 1/2 
and its z-component. By using the eigenvalues for J^, L^, and s^, we obtain I ■ s = 

+ 1) - + 1) - s{s + l)}/2 = -3/2 (1) for J = 3/2 (J = 5/2). As shown in 
Fig.l, therefore, E~ = — ^\ and = E^ + |A, where is the d-level energy 
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without SOI. In contrast to the rare-earth atoms, J = 3/2 and J = 5/2 states cannot 
be treated separately m 4d and 5d metal atoms smce the s-d hybridization potential 
is larger than the SOI. Therefore, we introduce the following SIAM for 4(i and 5d 
transition metal atoms with both p- and d-orbitals:^^) 

k,a kam fc,M,o=± 

+ E {^L^lP^m + h.C.} + Yl {VkAlAjM + h.C.} + ^ E 
kam k,(7,M,a=± M=iM' 

(2-1) 

Here, c^^ is the creation operator of a conduction electron with spin a = ±1. [d^jY 
is the creation operator of a d-electron with total angular momentum J = 2ibl/2 and 
z-component M {-J < M < J). 

Perm is the creation operator of a p-electron with 
angular momentum m(— 1 < m < 1). Sk = k"^ /2m is the energy for the conduction 
electrons, and (E^) is the localized d (p)-level energy. Here, we omit the atomic 
SOI for p-electrons since it is much smaller than \fi — E^\. V^j^j^ and are the 
mixing potentials, which are given by 

VkMa = ^y<i E ^'nl^Yrik), (2-2) 

VL = ^ypYrik), (2-3) 

where a*f^ is the Clebsch-Gordan (C-G) coefficient for J = 2 ± 1/2 and Yf^{k) is 
the spherical harmonic function. We will show that the phase factor in 1^'"(A;) and 
are indispensable to realize the SHE. Here, the C-G coefficient is given by 



o: 



mcr 



a{{5/2-Ma)/5y^U^^M-a/2 for J = 3/2, 



j^^+ = {(5/2 + Ma)/5}^/2^„,M-a/2 for J = 5/2. (2-4) 

In the present study, we neglect the crystalline electric field for each of the J = 3/2 
and J = 5/2 states to avoid too complicated expressions. We put h= 1 hereafter. 



L= 2, S=l/2/ 



J =5/2 



E- 



y=3/2 



spin-orbit 
interaction 

Fig. 1. The atomic SOI \l ■ s leads to a splitting between d-level with total angular momentum 
J = 3/2 and J = 5/2, whose energies are given by E' = E° - 3A/2 and E+ = E° + X, 
respectively. Here, E^ is the d-electron energy without SOI. 
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To discuss the scattering problem, it is useful to derive an effective Hamilto- 
nian for the conduction electrons by integrating out the d and p electrons in Eq. 
(^2-1). 36), 37) rj-ij^g obtained Hamiltonian is given by 



where 



k,k',(7 
Jka,k'a'^^ka^k'a' ^ 



"^kM 



J, 



ka,k'a' 



+ 



ka 

E 

k ,k' ,a ,ct' ,a=± 



— i — [Vw T 

jj^p / J km^ k'm' 

^ m 



^+ ^ ^kMa i^k^Ma') 

47rJ± Yl at{^a^^^^^,Y^\k) [y^^' {k') 

Mmm' 



(2-5) 



(2-6) 



(2-7) 



Jp = \Vp\^/{^l - EP) and J± = \Vd\'^/{^l - E^). Here, E^ = E^ + Rei7'^; S'^ is the 
d-electron self-energy due to the Coulomb interaction V^. According to the scaling 



-< 



A;, cr 



7± 



/ M \ U 

0---0 



' a' 



Fig. 2. Diagrammatic expression for J. 



theory, E^ approaches the Fermi level as the temperature decreases due to the 
Kondo effect: | J-|-| is strongly enhanced near the Kondo temperature Tk, and below 
Tk, J±N{0) ^ 1 due to strong resonant scattering, where A^(0) = mkp/2'K'^ is the 
density of state of the conduction band per spin. We also assume that JpN {{)){= 
-tan(5i/7r) < 1 since \^ - EP\ ~ 0(leV).35) 

Here, we explain the spin-orbit polarization due to the SOI, which plays a sig- 
nificant role in the intrinsic SHE.^^^ Using the relation (Z • = 1 (—3/2), the 
spin-orbit polarization ratio of d-electrons at the Fermi level is given by 
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N±{0) are the DOS for J = 2±l/2 state, which are given by N+ = EmL-5/2 Img^ (0)/7r = 
f Im^f (0) = miO)^J^, and iV_ = j:f_.^^^lmg^{0)/7r = |lm5^(0) = 4N{0)^J^, 
respectively. Here, g±{0) is the local Green function for d-electrons, which is given 

by 5^(0) = ^Ek^^^^GtiO)^ = ^Jy^^Z^Gt{0)5M,M'- From the 
above expression, note that {I ■ s)^ vanishes when = J^. 

§3. T-matrix and Current Vertex Correction 

In this section, we study the scattering problem. In both models given by Eqs. 
(2-1) and (2-5), the T-matrix due to the c-d resonant scattering is equivalent; it is 
given by 

'^ka,k'a' = '^kcr,k'a' + ■^kcr,k'a' 

TV X] (-^fca.fciai + '^fc<7,fciai)^fcl'^fcl(7i,fc'(T'' (^'l) 

fcl,CTl 

where its diagrammatic expression is shown in Fig. 3 (a), N is the number of fc-points, 
and G'^(£) = (e + — £k)~^- In this T-matrix, the terms containing both J^^y.^^^ 
and Jj^xax k'a' shown in Fig. 3 (b) vanish identically after ki and ci-summations as 
follows: 

fcl,cri o"i 

= Sm,M' XI ^ma]_(^mai = 0- (3-2) 



0-1 



Moreover, the term containing both J^„j^i„i and J^y given in Fig. 3 (c) vanishes 
identically due to the orthogonality of spherical harmonic functions. Then, the 
solution of Eq. (3-1) for k = k' and cr = o"' is simply given by 

rrid I \ 2[/ — 3T_^ , . 

T'ka,ka\^) = ^ TT + 7 TTT' (^'^j 

1 - J^g{£) 1 - 

where we have used the relations Y.M \^kMa^ = "A^d!? (SlV'dl^), and J2ka ^kMa 
X Glie) {V^M'aT = \yd\'^9{£)^MM'- oi^) = T,k ^ki^) IS the local Green function. 
Assuming approximate particle-hole symmetry near fj,, we put g^{0) = g{+i6) = 
—iTTN{0). Then, the quasiparticle damping rate in the T-matrix approximation is 
given by 

7d = -'^impImT^^fc^(O) 

7rA^(0)j2 ^ 7riV(0)4 



^''''"Pl + (7r7V(0)J_)2 + 3^i-Pi + (7riV(0)J+)2 
7_+7+, (3-4) 
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(a) 



= < X ^ + <X < X < + -^^r^ X ^<r^ 

J- J- J- J- J- J- 



+ < X — ^ 



< X < X < 



-^r^ X >(r^ 

J+ J+ J+ 



(b) 



(c) 



k, ctX ^ X*:': ct' 








7+ 



, '^l fc2 

— ^ 



kj.cr' 

k,ay4r- 



T± 
•-'ka.k^a' 







^k',a 



Fig. 3. (a) Diagrammatic expression for the T-matrix due to c-d resonant scattering, (b) and 
(c) Diagrams that vanish after fci -summation, (d) Expression of the Ward-Takahashi identity 
connecting between the current vertex correction (CVC) and the self- energy Sk- (e) Type I 
of CVC is obtained by differentiating the Green function, and (f) type II CVC obtained by 
differentiating the impurity potential. The former (latter) process gives the skew scattering 
term (side jump term). 



where njmp is the impurity concentration. Note that Eq. (3-4) is exact if riunp ^ 1- 
Here, we discuss the current vertex correction (CVC) due to a single impurity, 
which is the origin of the extrinsic terms. The total current in the presence of 
the CVC is given by the Ward-Takahashi identity as follows, which is shown in Fig. 3 
(d):38) 

Jk = Vk + (3-5) 

where the second term in the right-hand side is the CVC, which is required to satisfy 
the conservation law. Ej, = n^m^Tj^^h is a self-energy in the T-matrix approximation. 
As shown in Fig. 3(e) and (f), two types of CVCs are derived from the Ward identity. 
The type I of CVC is obtained by differentiating the Green function as shown in 
Fig.3(e), and the type II of CVC is obtained by differentiating the impurity potential 
as shown in Fig. 3(f). We will show that the skew scattering term arises from the 
former type of the CVC and the side jump term arises from the latter type. In 
contrast, the CVC due to the local impurity potential is negligible in the intrinsic 
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§4. Skew Scattering Term 

In this section, we study the skew scattering term using the Unear response 
theory. Initially, we consider the cases A^(0)|J±| ^ 1 and \Jp\ ^ |^±|) where the 
Born approximation is valid. In analogy to Refs. 21) and 35), the lowest skew 
scattering term is given by the following second Born approximation, 

ss(2ndBom) _ 6 1 \^ OSk' 2 

k,k'cr ^ y 

x|G^'(0)|^{TSt(0)J^,,, + c.c.}, (4-1) 

where — e (e > 0) is the electron charge. Its diagrammatic expression is shown in 
Fig.4 (a). Here, 




h' a' .h<y 



Fig. 4. Diagrammatic expressions for SHC induced by skew scattering (a) within the lowest order 
(extended Born approximation) contribution, and (b) in the T-matrix approximation with fuU 
order diagrams. 



fcl,(Ti,Q:=± 

= ^''(O)^ JL,fc'<x, (4-2) 

is the second-order term of the T-matrix; the first-order term in J^^^i^ does not 
contribute to <Tg|j up to the first-order term in J±?'^^ Note that any diagram that 
contains the part shown in Fig. 3 (b) vanishes identically. The retarded Green func- 
tion is given by G^(0) = (/i — + 27)""*^, where 7 represents the quasiparticle 
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damping rate. Here, we put 

7 = 7d + 7o, (4-3) 

where 70 is the damping rate due to nonmagnetic scattering, such as c-p scattering 
(7p = 37rninipA^(0) Jp) and the scattering due to disorders. The charge current is 
given by jf' = —edsk/dk^, where // = x,y. The spin current is then given by 



jf, = {a/2)deu/dk 



At- 



First, we consider the angular integration in Eq. (4-1), which is given by 



■'Ll,'A'.l,S(k,k')) (44) 



where f(k,k') = — ^— = sin 0^ cos 0^ sin 0^/ sin </)fc/ , and {■ ■ ■)o denotes the 
okx oky 

average over the Fermi level, which is defined as (A(k, k')\ = [ — f ^ A(k, k'). 

\ / Q J (47r)^ 

Since /(fc, k') = 2t: - l^f ^fc)} {>"f ^fc') + Yiik')] /3i, angular integration 

such as / dOkY-i^ '^^'^{k)Yi'^{k)Y-[^^{k) appears in Eq. (4-4). This integral is finite 
only when / = 1,3. Therefore, the interference of the d {I = 2) and p [l = \) partial 
waves is essential for skew scattering. After performing the angular integrations, 
Eq. (4-4) is given by 

E f (4V.4^, k'))^ = -'-J,{J^ - J„). (4-5) 

iT,a=± 

Using the relations |G^(0)|2 ^ ^6{^i - e^) for small 7, ^Zk, G^^i^) = g'^iO) = 
—inNiO), and the spin-orbit polarization {I ■ s)^ given in Eqs. (2-8) and (4-1) is 
transformed into 

^ss(2ndBorn) ^ ^ ^..^^ J^iil^k (3^ + 2j!)iV(0) 

- 1 (4.6) 



N"^ ^ dk dk 



7 



k,k' 

Since 7_ = 27rnimpiV(0) (7+ = 37rnimpiV(0) for J = 3/2 (J = 5/2) in the Born 
approximation, (Tgy is given by 

s.(2ndBorn) ^ ^J_j u4 ^i . 7d_ (4.7) 

27r 307r2 ^ ""^^ ^ ' 

From the above expression, Uglj vanishes when J+ = J_. 

Now, we derive the skew scattering term using the T-matrix approximation, 
which gives the exact result for nimp ^ 1- In this case, ^"^^^^^^(0) in Eq. (4-1) 

is replaced with the full T-matrix T^^^y,{Q) = Za=±Ti^alZ(0) {l - g^{0)Ja,y\ 
where the first-order term in J^^ has been dropped. The diagrammatic expression 
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for cjglj is shown in Fig. 4 (b). The angular integration in the T-matrix approximation 
can be performed as 



Jrt 



7riV(0)j2 



7rA^(0)j2 



1 + (^iV(O) J4 



1 + (7riV(0)J_)' 



(4-8) 



Then, the nonperturbative expression for (Tg|j with respect to J± obtained by the 
T-matrix approximation is given by 



1 



2tt 307r2 



X n 



imp 



37riV(0)4 +27r7V(0)j2 
[1 + (^iV(O) J+)2] [l + (^iV(0)J_)2] 



(4-9) 



where (Z • s) ^ is given in Eq.(2-8). In the case of A^(0) J+(_) » 1 and A'"(0) J_(+) ^ 1, 
the quasiparticle damping rate in Eq. (3-4) is given as 7+(-) • In the case of 7^ = 7 
and /cp = Tr/a, furthermore, the skew scattering term is simply given as 



27ra 10 
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(4-10) 



Here, n = /cp/37r2 is the density of the conduction electrons. Thus, dgy oc 7^/7^ and 
therefore, considerably large SHC can be realized due to the skew scattering term 
(Tgg in the low resistivity metals. 

§5. Side Jump Term I: for {I • s) ^ = -3/2 or 1 



In this section, we derive the side jump term for J = 5/2 and 3/2 indepen- 
dently under the assumption that the separation between these two states is complete 
(J+ ^ J- or J- ^ J+). However, when this separation is incomplete, the cross- 
term arises. We will discuss the contribution of the cross term in the next subsection. 
First, we consider the case where the Born approximation is valid: iV(0)|J±| ^ 1. 
In contrast to the skew scattering term, Jp is not necessary for the side jump term. 
The lowest order side jump term is given by the Born approximation. In the present 
model, it is given by a^^^°^"^ is given by 



sj(Born) 
'^SH 



(J = 5/2(3/2)) 



-n 



2tt 



imp 



^2 



k,k' 




+ Sm 



+{-) 



kMo 



dk. 




X |Gf (0)p 



E+(-) 



(5-1) 



whose diagrammatic expression is shown in Fig. 5(a). We stress that J~^^ ^ is neces- 
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Fig. 5. (a) Diagrammatic expression for the side jump term within the lowest order (Born approx- 
imation) contribution and (b) in the T-matrix approximation, (c) The diagram without J^'"^ 
vanishes identically. In contrast to the skew scattering term, Jp is not necessary for the side 
jump, (d) The diagram which vanishes identically unless Jp is taken into account. 



sary for finite SHC since the term without J"'"^"^ vanishes identically as discussed in 
section 7. In the present model, the charge current operator is given by = —evkfi, 
where — e (e > 0) is the electron charge, and 

^'"/^ = gfc^^fc^^fe'^ + I dk^ VkM + i^-c.j . (5-2) 

Next, we explain the Sz-spin current operator j^. In the present model, Sz is given 
by 

= E f 4.^.. + E S^t\d-J^^)U-J^^ + J: S',,idlt')'dtt\ (5-3) 

o- M M 

where Sj^~^^ is the spin Sz for only the J = 3/2 (5/2) state given by SjJ^^ = 

r j^/j f_L^1 ^ 

J2ma^ '^rna , and S'^j is that of the cross-term for J = 3/2 and 5/2 states 
given as 5^,^ = J2ma ^^ma ^ma ■ Straightforward to show that Sj^/^' = — ^ 

1 /2 

(f ) for J = 3/2 (5/2). Moreover, = {l - (^)^} for both J = 3/2 

and 5/2. We now discuss the case where J = 3/2 and 5/2 states are completely 
separated (i.e. J+ » J_ or J_ » J+). Then, S'j^^ in Eq. (5-3) is negligible, and the 

spin current = j'^fc^, Sz^ /2 is given by 

± 



Here, we study the velocity due to the c-d mixing potential VkMa given as^^^ 
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= v^, + vl (5-5) 

is the anomalous velocity given by the fc-derivative of the phase factor aM,a = 
exp {i (M — ^) 4>]^] in VkMa-^'^^'^^^ Since oc ky and thus v%{dek/dky) / 0, the 
anomalous velocity gives rise to the large SHE and AHE in heavy fermion systems. 
In contrast, oc k^ gives a normal velocity. Using the following relationships, 

ky Isin6'sin0 deu dsk . 

,2 I > 2 = 7 -2/1 ' -W = i^smgsmc/), (5-6) 

kj. + kfi k sm^ 6* afc^^ ok 

Eq. (5-1) is transformed as 

ai^-'^V = 5/2(3/2)) 

e 2 ^-^ ae 



-rii 



Y: ^ImG«(0)|G^(0)p 



fc,fc'(7,(T' 

2 



Vyl+^"V+^"Vy+^~h* I J+^"^ sin 

\ M I 



(5-7) 



where A'^}^'^ = 1{M - ^) + S^j- and k = \k\. First, we derive a^^ for 
J = 3/2. The angular integration in Eq. (5-7) is given by 

(^^MaykMaiyk'Ma'TJk'a'M ^^^^ '^^ ' ^^'^^ 

Note that this term includes the term calculated as J2a' {(^k' Ma'T'^k' M' a') n , ~ 
\V-\^5mm' ■ Using the following relations for J = 3/2, 



E^'lWl' = l^rfl'(l + 6sin2e), (5-9) 

Ma 

E^'lW = 4|^dP, (5-10) 

Ma 

E^^^l Wl' = 2|V^d|' (1 - 3sin2 0) , (5-11) 



Ma 

the angular integration in Eq. (5-8) is performed as 



E^M-lWl'sin^'^,} =-^j'-. (5-12) 



' \Ma I 



Then, using the relations |G^p ^(5(/i — efc) for small 7, Eq. (5-7) is transformed 
as 

1 / , N e 18 a 1 V 1 dsk vr^ ^, , 

'(J = 3/2) = Y^impJ% E aFT'^^'^ - '^^^^^ - 

1,1,' ^ 



^SH 
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By performing k and fe'-summations, figjj for J = 3/2 is given by 

^^T'^-XJ = 3/2) = -l^-imp^^^iV(O)^ (5-14) 

|^(Z-s)_. (5-15) 



2'Ka 5 7 

Here, {I ■ s)_ = —3/2 and we put kp = vr/a, where a is a lattice spacing. 

In a similar way to the calculation of cjgjj for J = 3/2, we derive cjgjj for J = 5/2. 
In this case, using the following relations, 

Y.^>kMf = ll^dWl + lOsin^e), (5-16) 

Ma 



E->>^Mf = m\', (5-17) 

Mo- 

T.Ma\V^,,f = 3|y,|2 (1 + 2sin2 0) , (5-18) 



Ma 

the angular integration is performed as 



^\Y1 "^MaKMai^k^Ma'T'^k'a'M si^^ 



\ M, 
a,a' 



Q 

= TSw-{T.^Ma\y^Ma?-^'^u) ='-^jl. (5-19) 

' XMa I Q 

Then, by using (Z • s)+ = 1 and 7+, dgy for J = 5/2 is given by 

c4^-^{J = 5/2) = ^ni„p^4iV(0)^ (5-20) 

^^^^-{l-s)+, (5-21) 



27ra 15 7 



where we put /cp = vr/a. 

Now, we derive the nonperturbative expression of the side jump term using 
the T-matrix approximation. In this case, J^i^^i^^ in Eq- (5-1) is replaced by the 
following full T-matrix: 

;+(-) 

rpdR+{-) _ "^kaMa' ,^ onx 

^Ha,U'a' -1_(^(0)J^(„))2- ^^"^^^ 

The diagrammatic expression for (Tg^ is given in Fig. 5 (b). In this case, the angular 
integration in eqs. (5-8) and (5-19) is replaced as 

T— ^ E {AMaV^MaiVyMa'TTt^ia' <t>k) 
' Ma 
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5 l + (^iV(0)J„)2 

11 4 

5 1 + (7riV(0)J+)2 



for J = 3/2, (5-23) 
for J = 5/2. (5-24) 



Thus, (Tgjj in the T-matrix approximation is given by (1 + ('7rA^(0) J-t)^) times eqs. 
(5-14) and (5-20). Therefore, the expression for the side jump term in eqs. (5T5) 
and (5-21) is vahd beyond the Born approximation, by considering 7-1-. 

§6. Side Jump term II: for general (Z • s)^ 

In the previous section, we have derived dgjj in the case where J = 3/2 and 
5/2 states are completely isolated. In this section, we discuss the case where these 
two states are not isolated, where the cross-term between these two states appears. 
Therefore, the contributions of these terms are indispensable for a more accurate 
calculation. Note that the cross-term does not exist for the skew scattering term 
cjglj because of the relationship in Eq. (3-2). Now, we derive the cross-term for the 

• 1 • , si-cross T ii T-> ■ i- si-cross(Bom) . . , 

Side jump term aJ^ . in the Born approximation, cjgjj is given by 




Fig. 6. Diagrammatic expression for the cross-term of the side jump term o-g^''°^^{a) within the 
lowest order (Born approximation) contribution, and (b) in the T-matrix approximation. 



sj-cross(Born) _ 6 1 f ^ f , Q' \ ^^kMu \ ^^k 

- 27r'''"^PiV2 2^ I2I2 V dk. ] dk, 

M,M',a=± 

X \G^m'-^ [(^fcW)*Gi^'(0)4V,fc. + {R^ A)] , (6-1) 
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where S'j^ is given in Eq. (5-3). Its diagrammatic expression is shown in Fig. 6(a). 
Here, the spin current operator is given by 

~ 2^ 2 dk 

a ^ 

+ + (6.2) 

We remind the readers that the side jump term originates from the anomalous ve- 
locity due to the c-d mixing potential. Using Eq. (5-6), Eq. (6-1) is transformed 

as 

sj-cross(Bom) 6 V~~^ 2 ds^ „R /„n 

k,k'a,a' , 



X 



\ M / 

X -^fc'a'.fcaSin^'/'fc. (6-3) 

where A'j^^ = ^(M — §)(^ + S'j^). As is the case with the previous subsection, the 
angular integration in Eq. (6-1) is performed as 

E ( E^A/'^^fcM.(^fcW)V,V,,.sin2</,A 

= E _\a^2 (jlAi.ykM.iysM.r ^>^) 



= -^(4-J!). (6-4) 
After taking the fc-summation in Eq. (6-3), ^^j^'^''"*'^^^™") given by 

Therefore, according to eqs. (5-14), (5-20) and (6-5), the final expression for the side 
jump term including the cross-term in the Born approximation is given as 



ir"' = i"<.,{y4-g^^}A'(0)^f 
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Now, we derive Cg'jj*^™'^'^ in the T-matrix approximation. In this case, J^i^ji jr,^. in 

Eq. (6-1) is replaced with the full T-matrix T^/^l^J- The diagrammatic expression 

for (Tgjj"°'^'^ is shown in Fig. 6. As is the case with the Born approximation, the 
angular integration in the T-matrix approximation can be performed as 



a=± ^ \ Ma I 

\ ( Jl Jl 



5 Vl + (^A^(0)J+)2 l + (7riV(0)J_)2 



(6-7) 



Therefore, the nonperturbative expression for the side jump term in the full T-matrix 
approximation is given as 



e2 fcF 1 / 37riV(0) J| 27rA^(0) J2 

2^3^ *)a'V7''^°'p \ [1 + (^iV(O) J+)2] [1 + (7riV(0)J_ 



In the case of A^(0) J_|_(„) ^ 1 and A^(0) J_(_|_) ^ 1, the quasiparticle damping rate 
in Eq. (3-4) is given as 7+(_)- In the case of 7^ = 7 and /cp = 7r/a, the side jump 
term is given as 



The relationship oc {I ■ s)^ is derived only when the cross-term is taken into 
account correctly. 

Here, we discuss the quasiparticle damping rate 7 dependence of af^ and a^^. 
From eqs. (4-9) and (6-8), af^ oc 1/7^, a''^^ oc (7^)° for 7^ > 70. When a^x oc 1/7^, 
the spin Hall angle due to the skew scattering mechanism is independent of 7^, and 
therefore, independent of the impurity concentration. This is consistent with the 
recent experiment by Niimi et al.^'^^ 

Finally, we discuss the ratio between (Tg|j and fjgjj. By using Eqs. (4-9) and 
(6-8), the ratio <Tgjj/<Tg|j is given by |f where n = /cp/Svr^ is the density of 

the conduction electrons. Therefore, Ugjj in cjsh exceeds (Tg|j in samples with high 
resistivity. 

§7. Negligible terms for the side jump term 



Here, we discuss the terms for the side jump term, which vanishes identically 
or negligible. In section 5, we stress that J~^^~^ is necessary for a finite side jump 
term. This is because the term without J'^^^\ which is shown in Fig. 7 (a), vanishes 
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identically since Im ^_^^(_)_^.^ = 0. Furthermore, we also stress that the diagram 
shown in Fig. 7 (b) vanishes identically after fc'-summation unless Jp is taken into 
account. At last, the side jump term we have studied is given by the diagrams with 
v^^ in Eq. (5-5) and We have dropped the diagrams with v'^ and as shown 
in Fig. 7 (c), since all of them are proportional to Wimp, and therefore negligible 
for riimp ^ 1. In summary, the extrinsic SHE is given by the side jump term when 
Jp = 0. Then, the dominant term of the side jump term, which is of order 7'^ and 
n? , is exactly given by Figs. 5 and 6. 




Fig. 7. (a) The diagrammatic expression for the side jump term without J'^^~' which vanishes iden- 
tically, (b) The side jump term with both T^^j^,^, and T^'^p^,^,, which vanishes identically unless 
Jp is taken into account, (c) The side jump term with dxV and dyV , which are proportional to 
fiimp. Thus, they are negligible for nimp ^ 1. 



§8. Discussion 



8.1. The skew scattering mechanism in two-orbital model 

In the previous section §4, we studied the SHC by using the Green function 
method. In this section, we discuss the skew scattering mechanism based on the 
Boltzmann transport theory. For this purpose, we study a simplified two-orbital 
model with M = ±3/2, assuming a strong crystalline electric field. In this model, 

KMa « Vd {V4Yi'^{k)6M,3/2a + ^2"'' (^) ^-3/2. } ■ 

Since V^jj'} oc Y^^ik) [Y^'^''{k)] oc e^^'^'^fe [g-si'^^fe] approximately, the second-order 
term of T'^ is simply given as 

j2^^2ia{cl>k-<t>k') _ g-2io-((/)fe-</)j,/) 

In the Boltzmann transport theory, the spin Hall resistivity due to skew scattering 
is pglj oc cr (fcx^^'U^C^'^ — ^ ^''7)) 1 where w represents the scattering probabil- 
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ity, which is proportional to n 



imp 



rj.di2)R jp 2 



due to Fermi's golden rule in the 



present model. According to Ref. 35), skew scattering occurs when the scatter- 
ing probability includes an assymetric component w^^{k — )■ k') oc Ime*^'^'=~'^'='-' oc 
(fc X k')z- In the present model, w^^ arises from the interference of the d and p 

scattering channel, w^^ G (^fci^feV'^A- fe' + '^•'^•) • ^SiCi, w^^{ka — )■ k'a) oc {J\ — 

^ V "1 * 

J^) X Ime*'^('^'=~'^fc'^ in this model since y contains the term Y^^{k) Y^^{k') oc 
g±i(0fe-0feO_ jj^ summary, a conduction electron with a hybridizes with l^ = —2a 
state due to the strong SOI, and therefore, the spin-dependent skew scattering prob- 
ability w^^{ka — ;> k'a) oc {J^ — )ime*""'^'^'=~'^'=') arises from the interference of the d 
and p angular momenta. Note that the skew scattering contribution vanishes when 
= Thus, the origin of the SHE due to the skew scattering mechanism is well 
understood based on the simplified two-orbital model. 

8.2. Estimation of the spin Hall angle 

Here, we estimate the magnitude of the spin Hall angle due to the skew scattering 

al?TT 2e , . , ^ ^ , e^n 

term: tanasH = — ^- The longitudmal conductivity is given by axx = , 

axx " 2m7 

where n = kp/Sir'^ is the density of the conduction electrons. Then, the spin Hall 
angle is given by 

ianasH = —JpN{0)^{l-s)^. (8-1) 

5 7 

Since tan (5i = — 7rA^(0) Jp,^*^^ Eq. (8T) can be rewritten as tan asH = 

5 7 

if |Jp|iV(0) < 1. 

Next, we discuss the condition for a giant extrinsic SHE due to 4(i and 5d metal 
impurities. Prom the expressions for the skew scattering term in Eq. (4-9) and the 
side jump term in Eq. (6-8), the condition for the giant SHA is (Z • s)^ ~ 0(1). This 
seems difficult to be realized in d-electron systems since J = 5/2 and 3/2 states are 
energetically close, that is, J+ ~ J_. In Fe, since the orbital magnetic moment is 
small, the orbital degrees of freedom are quenched. Therefore, {I ■ s)^ of Fe impurity 
in Au seems to be small. ^^^'^'''^ In contrast, the condition (Z • s)^ ~ 0{1) is realized in 
Ce and Yb atoms due to the almost complete separation between J = 3zb 1/2 state, 
and as a result, the giant SHE is expected to emerge due to Ce and Yb impurities. 

8.3. The relationship between intrinsic and side jump terms 

In this section, we discuss the relationship between the intrinsic term cjgj^ and 
the side jump term o"gjj. Since these two terms originate from the anomalous ve- 
locity, a close relationship between these two terms could be expected. However, 
quantitative calculation that supports this idea has not been performed. To com- 
pare these two terms quantitatively, we compare the coefficient given as c"^*^^j^ = 
(^^sii^'^^ / ^^{l ■ where a represents the atomic distance. 

According to Eq. (5-21), the coefficient of the side jump term for J = 5/2 is 
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given by c^^ = |. On the other hand, we have recently studied the intrinsic SHE 
based on the orbitally degenerate periodic Anderson model (OD-PAM) in Appendix 
A.^^^ The obtained intrinsic term for J = 5/2 is given as 

4H = ^^(i-«)M for J = 5/2, (8-2) 

where (Z • s)^ = 1. Therefore, the coefficient c™* for the OD-PAM is given by 
pint = ii. As a result, the coefficients c™* and c^J obtained from OD-PAM and 
SIAM are comparable. Furthermore, in the previous study, we calculated (Tg^ and 
{I ■ s)^ in various 4(i and f>d transition metals based on a multiorbital tight-binding 
model. The coefficient c'°* for Pt, for example, is estimated as ~ 3. This is 
because ~ 10^ and (Z • s)^ ~ 0.33 in Pt. For the other transition metals, c™* 
is estimated to be of order 1. Therefore, these two terms not only have the same 
origin, but also have quantitatively comparable contributions. 

Our study indicates that it might be difficult to distinguish the intrinsic mech- 
anism from the side jump mechanism when the d-ions are strongly disordered. For 
a more quantitative study, the study of both mechanisms based on a realistic model 
is needed. This is an important future problem. Another future problem is to study 
the ac Hall conductivity. Previously, the present authors studied the intrinsic ac 
anomalous Hall effect (AHE).^^^ Therein, the intrinsic ac AHC shows a prominent 
deviation from the Drude-type behavior. This coincides with a recent experiment 
by Kim et al.^^^ On the other hand, although no studies have reported the extrinsic 
ac AHE or SHE, both the skew scattering term and side jump term may show a 
Drude-type behavior. Therefore, ac SHE measurements will be useful to distinguish 
between the SHE due to the intrinsic effect and that due to the extrinsic effect. This 
is an important future work. 

8.4. Comparison with the study based on first principle calculation 

Recently, both the intrinsic and extrinsic SHEs have been studied based on the 
KKR Muffin-tin approximation.^^) The calculated skew scattering term based on 
the Kubo formula agrees well with that based on the Boltzmann formula. As for 
the side jump term, however, their obtained values seem to be much smaller than 
the side jump term based on the SIAM given in the present study. This difference 
seems to be due to a different definition of the side jump term. In Ref. 41), the type 
II CVC in Fig. 3(f) seems to be dropped. We expect that this is the reason why the 
side jump term obtained in their study is small in magnitude. 

8.5. Summary 

In summary, we have studied the extrinsic SHE due to transition metal impuri- 
ties based on J = 2 it 1/2 single-impurity Anderson model (SIAM). The analytical 
expressions for both the skew scattering term (Tg|j and side jump term have been 
derived. As is the case with the intrinsic term, it was found that both the skew scat- 
tering and side jump terms are in proportion to (Z • s)^. Therefore, (Z • s)^ ~ 0(1) 
is a necessary condition for giant intrinsic and extrinsic SHE. Moreover, we found 
a nontrivial close relationship between the intrinsic and side jump terms due to the 
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fact that these two terms originate from the anomalous velocity. We have shown that 
Ci^^t = a'^^/{l ■ s)^ in the OD-PAM is very close to C^j = a'^^/{l ■ s)^ due to dilute 
impurities described in SIAM. Furthermore, according to the previous study/^^ c™* 
in the 4d5d transition metals is estimated to be of the same order. fTo distinguish 
the intrinsic mechanism from the side jump mechanism, the ac Hall effect would be 
useful.'^^) 
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Appendix A 

Derivation of the intrisic term for J = 3/2 and J = 5/2 

In the main text, we have discussed the close relationship between the intrinsic 
term (Tg^ and side jump term o"gjj. Here, we derive the intrinsic term for J = 3/2 
and 5/2 based on the orbitally degenerate periodic Anderson model (OD-PAM). 

In the presence of the strong atomic SOI, the J = 5/2 level is about 2000 K 
higher than the J = 3/2 level. Therefore, we consider J = 5/2 and 3/2 states 
separately. We note that I ■ s = ^ + 1) — l{l + 1) — s{s + 1)] is given as 

l-s = l for J = 5/2, 
l s = -3/2 for J = 3/2. (A-1) 

Here, we introduce the following OD-PAM Hamiltonian: 

ka kM 

+ X] i^kMad^M^ka + VfcA/crcL^fcA/) 
Mka 

+ U <M<M'^ (A-2) 

where Cj[.^ is the creation operator of a conduction electron with spin o" = ±1. dj^^j 
is an operator of a d electron with total angular momentum J = 5/2 or 3/2 and 
z-component M {—J < M < J), efc is the energy for c electrons, E'^ is the localized 
d-level energy, and U is the Coulomb interaction for d electrons. VkMa is the mixing 
potential between c and d electrons, which is given by 

VkMa = ^ff^^^d Y ^^.YnOk, M, (A-3) 

m 
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where a!^^ is the Clebsch-Gordan (C-G) coefficient given by Eq. (2-4) in the main 



text and Y, 



is the spherical harmonic function. 



Here, the conduction and d-electron Green's function for the OD-PAM in the 
absence of the magnetic field are given as^''^'^^^ 



I VfcA/o-l 



G 



kMM' 



{uj)5mm' 



riOd 



(A-4) 



We note that G'j^^-{uj) = and G^j.j is the d-electron Green's function without 
hybridization given as 



riOd 
^kM 



1 



(A-5) 



Now, we consider the quasiparticle damping rate r{uj), which is mainly given 
by the imaginary part of the d-electron self-energy 17 (w). In the present study, we 
assume that F is diagonal with respect to M Fmm' = I^mm' and is independent of 
the momentum. Here, we perform a calculation of the SHC using this constant 7 
approximation. Then, the retarded (advanced) Green's function is given by 



UJ + IJ.- Ek 



\Vd\ 



oj + fi - E"" + 



uj + fi - E^ + {-)i-f 
-1 



(A-6) 



We calculate the intrinsic SHC UgJ^ based on linear response theory. The SHC 
at T = in the absence of the current vertex correction is given by CgJ^ = crgjj + a^-^ , 
where 



^11 



— VTr r^S^RJC/VA 
2ttN ^ V 



i.j=0 



(A-7) 



47riV 



Here, Cgjj and Ug^jj represent the Fermi surface and Fermi sea term, respectively. 

In the present model, the charge current operator is given by j'^ = —evk/j,, where 
— e is the electron charge and 



dVkMa t J ,1, 
Qf^ VkM + h.c. 



(A- 
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Next, we explain the s^-spin current operator j^. In the present model, Sz is given 
by 

^z = Y^ ^cl^Cka + SMd\j^dMk, (A- 10) 

o- M 

where Sm = Ylma f ['^mo-]^- It is straightforward to show that Sm = (x) ^^'^ 
J = 3/2 (J = 5/2). Then, the spin current is given by 

-.S _ \ ^ Q- (jgfc J 



2 3^:^'^-'''^" 



SHE originates from the anomalous velocity given by the fc-derivative of the phase 
factor in VkMa- We note that the terms composed only of d^SkdySk vanish identi- 
cally. 

Here, we calculate the SHC by neglecting CVC according to Eq. (A- 7), using 
eqs. (A-9) and (A-11). and are composed of the conduction electron term 
dsk/dk^ = df^Ek. and the hybridization term d^V^. As studied in Ref. 17), the 
relationship cjgjj ( i.e. Ugj^ » '^sh) obtained for small 7. 

According to Eqs. (A-7), (A-9) and (A-11), the Fermi surface term due to the 
product of dfj,V and d^Sk is given by is given by 



a, 



kMa 



^^^.MJGf(0)PGr(0)+c.c.' 



(A-12) 



First, we confine ourselves to the case for J = 5/2 state. Then, by using the following 
relationships for J = 5/2 

Y.M'\VkMa\' = ^ (1 + lesin^e) , (A.13) 

Ma 

Ya^\VkMa\^ = 2\Vd\^ (A-14) 
Ma 

^MalVkMal^ = (l + 2sin2^) , (A- 15) 

Ma 

ky _l singsin0 ^_9£fc 
kl + kl' k sinH ' dky - dk ^^^^^^"'P' 



Eq. (A-12) is transformed as 
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where k = \k\. 

Here, we analyze Eq. (A- 17) when 7 is smah enough. In this case, 

(0) =L-eu- . ' ^ 

i^i - Ek + ir^r\ (A-18) 

where efc = £fc H 7^7, and Fc = , 7. Since 7/(x^ +7^) = t^S(x) for smah 

fi — E"' [fi — E''-Y 

7, we obtain the fohowing relationship: 

|Gf(0)|2 ^ 



{^,-EkY + r^ 

(A.19) 

Substituting the above equation into Eq. (A- 17), we obtain the following relationship 
for small 7: 

J e 227r 1 dsk 



0"SH 



27rA^ 15 ^ k dk 

k 



Now, we approximate the conduction electron as a free electron. Then, cxgy for 
J = 5/2 is given by 

I 11 kp e 11 

^SH = ^15 2^'^^^=2^15'^"^ 

^^:Nfs{1-s)^, (A.21) 



27ra 15 

where (Z • s)^ = 1, a is the lattice spacing and A'^ps represents the number of large 
Fermi surfaces. The second line in Eq. (A-21) is obtained by putting kp = n/a. 

At last, we derive the intrinsic term for J = 3/2. In this case, we use the 
following relationships: 

Y.M'\VkM.\' = ^ (1 + 6sin2 9) , (A.22) 

Ma 



5^cT'|^fcAf.|' = 2|yd|2, (A-23) 

Ma 

J2McT\VkMa\^ = \Vdf (1 -Ssin^^) . (A-24) 

Ma 

Then, we can perform the calculation of CgJ^ in a similar way to J = 5/2. As a 
result, cjg^ for J = 3/2 takes a negative value as 

T 3 fcir ,^ e 3 

= -^5 2^"^^^ = -2^5^^^^ 

lNFs{l-s)t,. (A-25) 



27ra 5 

Here, {I • s)^ = —3/2, and we put kp = vr/a 
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Appendix B 

Accurate expression for the side jump term based on a single-impurity 

Anderson model for rare-earth atoms 

In our previous study, we derived the side jump term based on the single- 
impurity Anderson model. Therein, the definition of the spin current operator 
was made by mistake. Instead of defining 



it should be defined by using the spin operator s^ as is the case with Eq. (5-3) in 
the main text. Then, the spin current operator is defined as 

-■s _ \ ^ ^ d^k t 

o-M I ^^ ) 

where 5m = f [omo-]^- is straightforward to show that Sm = {+^) for 
J = 5/2 (J = 7/2). By using the following relations for J = 5/2, 

E ^VfcMj' = + 16sin2 0), (B.2) 

Ma 



E^'IWP = 6|V7P, (B-3) 

Ma 

E^^l = 3|V/P (1 - 4sin2 0) , (B.4) 



Mo 



and for J = 1/2 



E ^'l^fcMJ' = + 30sin2 0), (B-5) 

Ma 

E-V^J^ = 8|y/P, (B.6) 

Ma 

E^^^l^fcA/J' = 4|l^/l'(l + 3sin2e) , (B.7) 

Ma 

the final expression for (Tg^ is given as 



^SH 



e Qkp 
2117 IT 7 
e 3 

2^7 
e 15 kp 7/ 
211 14 TT 7 
_^5 
2^7 



(Z • s);, for J = 5/2, (B- 



(Z-s)^ for J = 7/2, (B- 
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where spin-orbit polarization {I ■ s)^ is given by (Z • s)^ = —2 (3/2) for J = 5/2 
(J = 7/2). We note that only the coefficient of the final expression is slightly- 
modified compared with that in Ref. 28). 
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